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Abstract 

A D-dimensional gravitational model with Gauss-Bonnet term is considered. 
When ansatz with diagonal cosmological type metrics is adopted, we find solutions 
with exponential dependence of scale factors (with respect to “synchronous-like” 
variable) which describe an exponential expansion of “our” 3-dimensional factor- 
space and obey the observational constraints on the temporal variation of effective 
gravitational constant G. Among them there are two exact solutions in dimensions 
D = 22, 28 with constant G and also an infinite series of solutions in dimensions 
D > 2690 with the variation of G obeying the observational data. 
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1 Introduction 


Here we deal with D-dimensional gravitational model with the Gauss-Bonnet term. The 
action reads 

S=f d D zy/\i\{a 1 R[g\ + a 2 C 2 [g}}, (1.1) 

J M 

where g = gMNdz M <g) dz N is the metric defined on the manifold M, dirnM = D, \g\ = 
|det(0 M Jv)| and 

n rj jdMNPQ a td tjMN , rj2 (i o\ 

*—2 ~ P-MNPQP — 4KmNP- + OC 

is the standard Gauss-Bonnet term. Here an and a 2 are non-zero constants. 

Earlier the appearance of the Gauss-Bonnet term was motivated by string theory 
[D 121 El H E]. At present, the (so-called) Einstein-Gauss-Bonnet (EGB) gravitational 
model and its modifications are intensively used in cosmology, see [Q !7j (for D — 4), 
(siEuinuminunuiiuisuiT! and references therein, e.g. for explanation of accelerating 
expansion of the Universe following from supernovae (type la) observational data [181 USs 
m\. Certain exact solutions in multidimesional EGB cosmology were obtained in 0 -na 
and some other papers. 

Here we are dealing with the cosmological type solutions with diagonal metrics (of 
Bianchi-I-like type) governed by n scale factors depending upon one variable, where n > 3. 
Moreover, we restrict ourselves by the solutions with exponential dependence of scale 
factors (with respect to “synchronous-like” variable r) 

di(r) ~ exp (h l r), (1.3) 

i = 1,..., n; D = n + 1, with the aim to find solutions describing an exponential isotropic 
expansion of 3-dimensional flat factor-space, i.e. with 

h l = h 2 = h 3 = H > 0, (1.4) 

and small enough variation of the effective gravitational constant G, which is proportional 
to the inverse volume scale factor of the internal space, i.e. 

n 

G ~ J^Mt)]” 1 ~ exp (-Int r), (1.5) 

2=4 

where here and in what follows we denote 

n 

Int = h\ (1.6) 

2=4 

see [2, [221 22 M, ES1J and references therein. We call h l = di/ai as “Hubble-like” 
parameter corresponding to i-tli subspace. 
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In cosmological model under consideration with anisotropic “internal space”, we get 
for dimensionless parameter of temporal variation of G the following relation from (1.5) 
and (|1.6[) 


G 

GH 


= — Int/H , 


(1.7) 


where H is the Hubble parameter. 

As for the experimental data, the variation of the gravitational constant is allowed at 
the level of 1CT 13 per year and less. We use the following constraint on the magnitude of 
the dimensionless variation of the gravitational constant 


— 0, 65 • 1CT 3 < —— < 1,12 • 1CT 3 , (1.8) 

Gii 

which comes from the most stringent limitation on G-dot obtained by the set of ephemerides 

m 


G/G= (0.16 ±0.6) ■ lir 13 year- 1 (1.9) 

allowed at 95% confidence (2-cr) level and the present value of the Hubble parameter pTj 
(which characterizes the rate of expansion of the observable Universe) 

H 0 = (67, 80 ± 1, 54) km/s Mpc~ l = (6.929 ± 0,157) • HT 11 year~\ (1.10) 


with 95% confidence level. It shold be noted that the original result for Ho in [27] (for the 
Planck best-fit cosmology including external data set) was presented at 68% confidence 


(1-cr) level. In restriction (1.8) we use the lower allowed value for H 0 in (1.10) in order to 
obtain the confidence level more than 95% . 


Thus, we are seeking here the cosmological solutions which obey the relations (1.3)- 


(1.8), listed above. 

The paper is organized as follows. In Section 2 the equations of motion for D- 
dimensional EGB model are considered. For diagonal cosmological type metrics the 
equations of motion are equivalent to a set of Lagrange equations corresponding to certain 
“effective” Lagrangian [29] [30] (see also [9]l28)). In Section 3 some cosmological solutions 
with exponential behavior of scale factors satisfying the restriction (1.8) are obtained for 
two isotropic factor spaces and positive value of a = a^/oq. 


2 The cosmological type model and its effective La¬ 
grangian 

2.1 The set-up 

Here we consider the manifold 


M = M* x Mi x ... x M n , 


( 2 . 1 ) 
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with the metric 


g = me 2 ^ u >du ® du + ^ e 2 ’’* <“>„.*/ ® dy\ 


( 2 . 2 ) 


2=1 


where w = ± 1 , £* = ± 1 , % = l,...,n, and are one dimensional manifolds 

(either M or S 1 ). Here and in what follows M* = (ii-,u + ) is an open subset in M. The 
functions 7 (u) and /3 l (u), i = 1,... , n, are smooth on M* = («_,«+). 

For w = —1, £1 — ... — e n — 1 the metric (2.2) is a cosmological one, while for w — 1, 
£1 = — 1 , £2 — ... = £ n — 1 it describes certain static configurations. 

For physical applications we put M\ = M 2 = M 3 = M, while M 4 ,...,M n will be 
considered to be compact ones (i.e. coinciding with S' 1 ). 

The integrand in (1.1[), when the metric (2.2) is substituted, reads as follows 


where 


and 


^/\g\{oLiR[g\ + a 2 C 2 [g]} — L + -J-, 

(2.3) 

L = a\Li + a 2 L 2 , 

(2.4) 

Li = (-w)e-^ 0 GijFF, 

(2.5) 

L 2 = ~^e- 3 r '+*>G m P' 0 3 $ k $ l , 

( 2 . 6 ) 

n 

7o = £/?‘ 

2=1 

(2.7) 

$ij I? 

( 2 . 8 ) 

— l)(<5jfc — l)(^i/ — 1 — 1 )(5ji — 1 )( 8 m ~ 1 ) 

(2.9) 


are respectively the components of two metrics on M n [291130] . The first one is the well- 
known “minisupermetric” - 2-metric of pseudo-Euclidean signature: < v\,v 2 >= Gtjv\v 2 
and the second one is the Finslerian 4-metric: < v\, v 2 , U 3 , V 4 >= Gijkiv\v 2 v^vv s = 
(v*) G M n , where < .,. > and are respectively 2- and 4-lin ear symmetric forms 

on MA. Here we denote A = dA/du etc. The function f(u) in ( |2.3| ) is irrelevant for our 
consideration (see [291 ED]). 


The derivation of (2.4)-(2.6) is based on the following identities 

n n 

Gi,vV = J >‘) 2 ~ (£ 


,i \2 


( 2 . 10 ) 


2=1 


2=1 


G m vVv k v l = (]>>‘) 4 - 6(J>‘) 2 £( 




2=1 


2=1 


3 = 1 


+3(£( 


,i\2\2 


2=1 


+s(E3 i,< ) - 6 EZf 1 ’’)' 1 


( 2 , 11 ) 


2=1 


2=1 
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( 2 . 12 ) 


It follows immediately from the definitions (12.81) and (2.9) that 


Gij-uV = -2 J 2 

i<j 


V l V 3 , 


G m v^v k v l = 24 ^ 

i<j<k<l 


v l v 3 v k v l . 


(2.13) 


2.2 The equations of motion 


The equations of motion corresponding to the action (1.1) have the following form 


where 


&MN — ® 1 &MN ^2 &MN — 


^MN ~ RmN ~ V RgMNi 


£mn ~ 2 (RmpqsRn PQS — 2 R M pR 


■N 


1 


-2RmpnqR PQ + RRmn) ~ 2^ 2 9mn- 


(2.14) 


(2.15) 


(2.16) 


It was shown in j[3U] that the held equations (2.14) for the metric (2.2) are equivalent 
to the Lagrange equations corresponding to the Lagrangian L from (2.4). 


Thus, equations (2.14) read as follows 

waxGijfrfr + aie^Gijkifrp 3 i3 k j3 l = 0 , 

d[-2 wa 1 G ij e- l+ -'°P - /?*/)'] -1 = 0, 


(2.17) 

(2.18) 


i = 1,n. Due to (2.17) L = —w\e 7+70 a . 


2.3 Reduction to an autonomous system of first order differen¬ 
tial equations 

Now we put 7 = 0 and denote u = r, where r is a “sy nchro nous- like” variable. By 
introducing “Hubble-like” variables h 1 = f3 l , the eqs. (2.17) and (2.18) may be rewritten 
as follows 


wa\Gijh l h 3 + a2Gijkih t h 3 h k h l = 0 , 


—2 waiGijh 3 — -ot2Gijkih 3 h k h l 


J2 h ' 


i =1 


d 

dr 


-2waiGijh 3 - -a. 2 G ijk ih 3 h k h l 


-L = 0, 


(2.19) 


( 2 . 20 ) 
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i — 1,..., n, where 


( 2 . 21 ) 


L = —w(y.]G l 3 h l h 3 -a 2 Gij k ih l h 3 h k h . 

3 


Due to (2.19) L = wa\G i 3 h l h 3 . 


Thus, we are led to the autonomous system of the first-order differential equations on 
h 1 (r), h n {r) HUSO]. 


Here and in what follows we use relations (2.10), (2.11) and the following formulas 


G ij v i =v i -S 1 , (2.22) 
G ijkl v j v k v l = Sl + 2S 3 - 3SiS 2 + 3(^ 2 - S*y + 6Si(r/) 2 - 6(r/) 3 , (2.23) 

i = l,...,n, where S k = S k (v) = £" = i («*)*. 

2.4 Solutions with constant h l 


In this paper we deal with the following solutions to equations (2.19) and (2.20) 

h\r) = v\ (2.24) 

with constant v l , which correspond to the solutions 

P = v l r + fa, (2.25) 


where p l Q are constants, i — 1,..., n. 


In this case we obtain the metric (2.2) with the exponential dependence of scale factors 

n 

g = wdr <g) dr + E l Ble lv '" T dy 1 <S> dy\ (2.26) 

i =1 

where w = ±1, £; = ±1 and B, t > 0 are arbitrary constants. 

For the fixed point v — (u*) we have the set of polynomial equations 


GijV l v j — aG i j k iv l v 3 v k v l = 0, 


2 GijV 3 - -aG ljkJ v 3 v k v l ^- 


2—1 


^GijV l v 3 = 0 , 


(2.27) 

(2.28) 


i = 1,... ,n, where a = a 2 {—w)/a\. For n > 3 this is a set of forth-order polynomial 
equations. 

The trivial solution v = (v 1 ) = (0,..., 0) corresponds to a flat metric g. 

For any non-trivial solution v we have \ v% 7^ 0 (otherwise one gets from (2.28) 
GijV l v 3 = Er=i(^) 2 — (Z jLi V ' ) 2 = 0 an d hence v — (0,..., 0)). 

The set of equations (2.27) and (2.28) has an isotropic solution v 1 = ... = v n = H, 
where 

n(n — 1 )H 2 + an(n — 1 )(n — 2 )(n — 3 )H i = 0. (2.29) 
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For n — 1: 77 is arbitrary and for n = 2,3: 77 = 0. 

When n > 3, the non-zero solution to eq. (2.29) exists only if a < 0 and in this case 


[291130] 

77 = ±^= 1 (2.30) 

V\ a \( n - 2 )(n - 3) 

In cosmological case (u> = —1) this solution takes place when a 2 /ai = a < 0. 

The isotropic solution for n > 3 gives rise to a very large value of G/G = (n — 3)77, 
which is forbidden by observational restrictions. 

It was shown in [29] ED] that there are no more than three different numbers among 


3 Examples of cosmological solutions obeying the re¬ 
striction on the variation of G 


In this section we consider some solutions to the set of equations (2.27), (2.28) of the 
following form v = (H,... ,H,h,... ,h), where 77 the “Hubble-like” parameter corre¬ 
sponding to m-dimensional isotropic subspace with m > 3 and h is the “Hubble-like” 
parameter corresponding to /-dimensional isotropic subspace, l > 2. 

These solutions should satisfy the following conditions: 


1. mandatory: 

(a) 77 and h are real numbers, 

(b) 77 > 0, h < 0; 

2. desirable: 


(a) Int — (rri — 3 )H + lh < 0; 

(b) -0, 65 • 1(T 3 <-§h = -((m - 3) + |) < 1,12 • 1(T 3 . 

The first inequality 77 > 0 in the mandatory condition is necessary for a description of 
accelerated expansion of 3-dimensional subspace, which may describe our Universe, while 
the second inequality h < 0 excludes an enormous (of Hubble’s order) variation G/G for 
h > 0 and m > 3. 

The first desirable condition means that the volume scale factor of the internal space 
V(t) = 7?exp(((m — 3)77 + /h)r), where B > 0 is constant, decreases over time. This 
condition is a sort of weak extension of a possible restriction for rn — 3 coming from 
the unobservability of the “internal space” for all t > Tq. It is also desirable since the 
negative value of the para meter Int is more probable due to more probable positive value 
of G/G = —Int, see (1.9). 

2 At the moment we were unable to find solutions with three different real “Hubble-like” parameters. 
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The second desirable condition may be also rewritten by using the parameter V ar = 
|^| = |(m-3) + t|: 

V ar < 1,12 • 10 -3 , for Int < 0; (3.1) 

Var < 0, 65 • 1(T 3 , for Int > 0. (3.2) 

Here we consider the simplest case, when internal spaces (apart the expansion factors) 
are flat. The consideration of curved internal spaces will drastically change the equations 
of motion and may break the existence of solutions with exponential dependence of scale 
factors. Anyway the inclusion into consideration of curved internal spaces may be worth 
but need a special treatment, which may be done in a separate work. 


3.1 The dependence of “Hubble-like” parameters on m and l 

The total dimension of the considered space is D = n + 1 = (in + l) + 1, where we have 
m dimensions expanding with the Hubble parameter H > 0 and l dimensions contracting 
with the “Hubble-like” parameter h < 0. 

According to this, we rewrite the set of polynomial equations (2.27), (2.28), using the 
formulas (2.22), (2.23), as follows: 


H 2 (m — m 2 ) + h 2 (l — l 2 ) — 2mlHh 
—a(H 4 m(m — l)(m — 2)(m — 3) + hH(l — 1)(Z — 2)(/ — 3) 

+4 H 3 hm(m — 1 )(m — 2)1 + 4 h 3 Hl(l — 1 )(l — 2 )m 

+ 6 H 2 h 2 m(m — l)l(l — 1)) = 0, (3.3) 


m(l - m)H 2 - (1/2)Z/i 2 (1 + 21) + 2ZHA((3/4) - m) 

—a(H 4 m(m — l)(m — 2 )(m — 3) + H 3 hl(m — l)(m — 2) (4m — 3) 

+3 H 2 h 2 l(m — 1)(2 Im — 21 — m) 

+Hh 3 l(l - l)(4Zm -31- 2m) + hH 2 (l -l)(l - 2)) = 0, (3.4) 


1(1 - l)h 2 - (l/2)mH 2 (l + 2m) + 2mHh((3/A) - l) 

—a(h 4 l(l - 1 )(l- 2)(l - 3) + h 3 Hm(l - 1 )(l - 2)(4 1 - 3) 

+3 h 2 H 2 m(l — 1)(2 Im — 2m — l) 

+hH 3 m(m — 1)(4 Im — 3m — 21) + H 4 m 2 (m — 1 )(m — 2)) = 0. (3.5) 

Here we put for simplicity a = ±1 but keep in mind that general ct-dependent solution 
has the following form 


H(a) = H\a\~ 1/2 , h(a) = h\a\~ 1/2 . (3.6) 

Due to these relations the parameter G/(GH) does not depend upon a and hence 
our simplification is a reasonable one. For any solution (H, h) with a = ±1 we can find a 






proper a which will be in agreement with the present value of the Hubble parameter H 0 
(see (|1.10|)) 

H\a \~ 1/2 = H 0 . (3.7) 

Our numerical analysis (based on Maplesoft Maple) show us that (generically) there 
are 11 solutions of these equations (for m > 3 and l > 3). 

I) The first to mention is, obviously, the zero solution H\ = h\ = 0. 

II) Two other solutions are isotropic ones: 

1. if a = 1, then H = h = ± J (f 3+2im+m L 5t _ 5m+6) - ■ We are led to pure imaginary 
isotropic solutions obeying (2.29). For example, when m = 9 and l — 6 we receive 

(a) H 2 = h 2 = \J ^ ■ i; 


(b) H 3 = h 3 = 2 y 39 ' b 

2. if a = — 1, then H — h — ±^/ ^ 2 + 2 lm+m 2 _ 5 i_ 5m+ ^ ■ We are led to isotropic solutions 
(2.30). When m — 9 and l = 6 the solutions are: 


(a) H 2 -h 2 - 2 ^; 

(b) H 3 = h 3 = -\J±. 

Ill) For a = ±1 the remaining eight solutions are roots of the following two equations 
which are given by Maple: 


P(H ; m, l) = 64(m — 2 ){m — l) 2 (m — 2 + l)(l 2 m + lm 2 — 2 1 2 

+2 Im — 2m 2 )(—3 + l + m ) 2 ■ H 8 
=Fl28 (m — 1) 2 (—3 + l + m){l 2 m{l + m) 2 — 21(1 + m) 3 
+ (/ + m)(8/ 2 + lm + 2m 2 ) — 10Z 2 + 4/m — 4m 2 ) • F7 6 
+16(m — l)(5/ 5 m + 10/ 4 m 2 + 5 Z 3 m 3 — 6/ 5 — 38/ 4 m 

—49/ 3 m 2 - 17Z 2 m 3 + 32/ 4 
+lll/ 3 m + 75 Z 2 m 2 + 14/?n 3 — 70/ 3 — 130/ 2 m — 14/?n 2 
—8m 3 + 68Z 2 + 4/m + 8m 2 ) ■ 77 4 
=Fl6/(m — 1)(Z 4 + Z 3 m — 3/ 3 — 5Z 2 m + 5Z 2 + 8Zm — 71 — 2m) ■ H 2 

+Z 5 - 3Z 3 - 2Z 2 = 0, (3.8) 


P(/i; Z, m) = 0. 


(3.9) 


The second equation is obtained from the first one just by swapping parameters m 
and Z and replacing H by h. The solutions to eqs. (3.8), ( 3.9[ ) should be substituted into 
eqs. (3.3), (3.4), (3.5), in order to find the solutions (H,h) under consideration. 
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The closed-form expression for the solution in general case (for any m and /) seems to 
be very bulky. So, we use Maplesoft Maple to find solutions for certain m and / and test 
some general features of these solutions 4 : 

1. for a — 1 in common case we have two pairs of real and quad of complex solutions 
which differ in signs. (See footnote 4.) For m — 9 and / = 6 we receive: 

(a) ff 4 « -0.2597 and h A « 0.2826; 

(b) ff 5 « 0.2597 and h 5 « -0.2826; 

(c) Hq ~ —0.1610 and h 6 ~ 0.4004; 

(d) H 7 « 0.1610 and h 7 « -0.4004; 

(e) 77 s « -0.0913 - 0.0464 • i and /i 8 « 0.1456 - 0.0449 • i; 

(f) #9 ~ -0.0913 + 0.0464 • i and /i 9 « 0.1456 + 0.0449 • i; 

(g) ifio « 0.0913 - 0.0464 • i and /ii 0 « -0.1456 - 0.0449 • i; 

(h) iiii « 0.0913 + 0.0464 • i and /in «■ -0.1456 + 0.0449 • i; 

2. for a = — 1 (m > 3 and l > 3) there are no (extra) real solutions. (See footnote 4.) 
In the case of m = 9 and / = 6 we obtain: 

(a) H 4 « —0.2597 • i and h 4 ~ 0.2826 • i; 

(b) ff 5 « 0.2597 • i and /i 5 « -0.2826 • i; 

(c) ii 6 ~ 0.1610 • i and /i 6 ~ —0.4004 • i; 

(d) if 7 ~ —0.1610 • i and h 7 ~ 0.4004 • i; 

(e) ff 8 w -0.0464 - 0.0913 • i and /i 8 w -0.0449 + 0.1456 • i; 

(f) ff 9 « -0.0464 + 0.0913 • i and /i 9 « -0.0449 - 0.1456 • i; 

(g) ifio « 0.0464 - 0.0913 • i and h 10 « 0.0449 + 0.1456 • i; 

(h) if n « 0.0464 + 0.0913 • i and /in « 0.0449 - 0.1456 • i. 

It can be seen that none of the solutions for a = — 1 satisfies our mandatory conditions 
written in the beginning of this section. 

As for a — 1, the solutions III. 1 .b and Ill.l.d are real and H > 0, h < 0. It can be 
verified that in these cases Int = (m — 3)if + Ih < 0 . 

Now we have to calculate the variation of the gravitational constant. For m — 9 and 
/ = 6 we get: 

I/ar 5 = | Aj _| 5 = |(m - 3) + « 0.535826; 

t/ar 7 = |^| 7 = |(m - 3) + « 8.914741. 

4 Here we are led to some common features of the solutions just by numerical calculations for a restricted 
range of numbers in and l. 
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The first variation is lower: Var 5 < Var 7, for m = 9 and l = 6. This inequality seems 
to take place for any m > 3 and l > 3. At the moment a rigorous proof of this fact 
is absent while certain numerical calculations support it. Anyway, here we will focus on 
the solution ( H 5l /i 5 ), which we consider as a more interesting (for our applications) then 
(H 7l h 7 ). Further we will write H and h instead of H 5 and h 5 in common case. 

We can plot behaviour of parameters H , h, Int and Var , for example, keeping fixed 
m = 8 and m — 10 and raising l from 5 to 100 by 5. See Figure 1 and Figure 2. 


H(l) 



h(i) 

0 20 
0 -- 

-0,05 - 

- 0 , 1 - 

-0,15- 

- 0 , 2 - 

-0,25 - 

-0,3- 

-0,35 - X 



40 60 80 100 

J 1 1 1 I 1 1 1 I 1 1 1 L 

,.********* 
* * * * 


♦ for m=8 
x for m=10 


Figure 1: Behavior of “Hubble-like” parameters H and h for fixed m = 8 and m — 10 
while l is changing. 


Int(l) 

0 20 


Var(l) 




Figure 2: Behavior of the internal space parameter Int and the variation of G parameter 
Var for fixed m = 8 and m = 10 while l is changing. 
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3.2 The limiting values of H , hi, Int and V ar for fixed m < 9 


When m < 9 the internal space parameter Int remains negative, that means that the first 
desirable condition is satisfied for any l. The variation of G parameter is monotonically 
decreasing with the increase of l. Moreover, we get finite limits for H and hi as l —> +oo. 
In this subsection we obtain these and other limits (for Int and V ar) for fixed m < 9. 

Now let us rewrite (3.8) and (3.9) for H and hi keeping only the terms with higher 
degrees of/: 


64 (m - l) 2 (m - 2 ) 2 / 5 ■ H 8 - 128 (m - 2)(m - 1 ) 2 / 5 • H 6 

+16(m - l)(5m - 6 )/ 5 • H A - 16(m - l)l 5 ■ H 2 + l 5 = 0, (3.10) 

64(m — 1)(/ • h ) 8 — 128(m 2 — 2 m + 2)(/ • h ) 6 
+ (80?n 3 - 272 m 2 + 224m - 128) (l ■ hf 
— 16 m{m 8 — 5 m 2 + 8 m — 2 )(l ■ h ) 2 

+m 5 — 3m 3 — 2 m 2 = 0. (3.11) 


Solving these equations we find the limiting values: 


, TT 1 2m — 2 + \j2m 2 — 2m 

Inn H = —\ ---, 

i^-oo 2 V rn 2 — 3m + 2 


lirn h-l = — { - 

i —^oc 2 Vm — 2 


(m^/2m(m — 1) + 2?n 2 


- 2 , 


'2 ((m + 2 )y / m(m — 1) + m 2 v / 2)(m — l) 2 


— 1 ) 


—4m + 4) 


(3.12) 


(3.13) 


lim Tar 

Z—> oo 


(m 


3 ) + 


linp^oo h • l 

lim^oo H 


(3.14) 


We fill the Table 1 by calculated values for 3 < m < 9: 
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m 

lim^oo H 

lim^oo l ■ h 

linh-kxi Int 

lim Var 

3 

0.9659258265 

-0.7630807575 

-0.7630807575 

0.7899993318 

4 

0.6738873385 

-1.093021916 

-0.4191345775 

0.621965355 

5 

0.5462858555 

-1.352249104 

-0.259677393 

0.475350754 

6 

0.4709825726 

-1.574449592 

-0.161501874 

0.342904141 

7 

0.4199717390 

-1.772664074 

-0.092777118 

0.220912766 

8 

0.3825276619 

-1.953613607 

-0.040975297 

0.107117214 

9 

0.3535533906 

-2.121320344 

0 

0 


Table 1: The limiting values of H, hi, Int and Var parameters as l —> oo. 


For 3 < m < 8 the limiting values of Var-parameter are too large. Since Var{l) 
exceed the limiting values Var(oo) for 2 < m < 9 (see Figure 2 for m = 8) the restriction 
(1.8) on variation of G is not satisfied for m = 3,..., 8 and we are led to unphysical results. 
This is why we consider in what follows just the cases m = 9,10,12,.... 


3.3 Infinite series of solutions for m = 9 


Now we consider the case m = 9. We get the following relations. 


For 1 = 2679: 

H « 0.3531582111, 
h « -0.0007910955039, 
Int ps -0.000395588, 
Var ps 0.001120145. 

For 1 = 2680: 

H ps 0.3531583594, 
h ps -0.0007908005933, 
Int ps, -0.000395434, 
Var ps 0.001119706. 


For 1 = 2681: 

H ps 0.3531585062, 
h « -0.0007905059028, 
Int ps -0.000395288, 
Var ps 0.001119293. 

For 1 = 2682: 

H ps 0.3531586532, 
h ?s -0.0007902114318, 
Int ps -0.000395141, 
Var ps 0.001118876. 


We do not present here exact analytical forms of these solutions in radicals which are 
bulky ones. For example, the relation for the parameter H, when l = 2680, contains (17 
times) the radical V32839778319264444823234828568184005. 

The numerical calculations for fixed m — 9 gives us an evidence of monotonically 
decreasing behaviour of the function Var(l) for l > 2680 5 as well as the asymptotical 
relation: Var(l ) ~ A/l, as l —> +oo, where A > 0. See Figure 3. 

Thus, for rri — 9 there is an infinite series of admissible cosmological solutions with 
l = 2680, 2681,..., which satisfy all the conditions imposed. Any of such solution describes 
accelerated expansion of the three-dimensional factor space with sufficiently small value 
of the variation of the effective gravitational constant G. This variation may be arbitrary 
small for a big enough value of l. 

5 A rigorous analytical proof of this fact may be a subject of a separate work. 
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Var(l) 



Vax(l) 



Figure 3: Behavior of the variation of G parameter Var for fixed m — 9 while l is changing 
from 1000 to 11000 by 500 (left) and from 2670 to 2690 by 1 (right). 


The infinite series of solutions for m = 9 and l = 2680, 2681,... starts from the (special) 
total dimension D = 2690. For D < 2690 and m — 9 the solutions do not obey restriction 
(1.8) on variation of G and hence are not of interest for our consideration. 


3.4 Some solutions for m > 9 with minimal War-parameter 

When m > 9 the internal space parameter Int becomes positive. As l can only be a 
natural number we should look for a value of l which gives the minimal magnitude of the 
variation of G parameter Var. Below we present the calculated values of our parameters 
(H, h, Int — (m — 3) ■ H + l ■ h and Var = |(m — 3) + ^|) for each of considered cases: 

1. for m = 10: 

the variation of G parameter is minimal for 
1 = 31, see Figure f, 

H « 0.2996055415, 

h « -0.06764217686, 

Int « 0.000331307, 

Var » 0.001105812. 

This case is not of particular interest. The 
radical forms of the solutions are too bulky, 
so we approximated them. Nevertheless the 
variation of G parameter is out of the al¬ 
lowed domain and the “internal space” , T1 . ,. r 

1 igure 4: 1 he variation of G parameter tor 
parameter is positive. m — 10 


Var(I) 
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2. for m = 11: 


the variation of G and the “internal space” 
parameters are zero for 1 = 16, see Figure 0,M “ 
5, 0 , 12 - 



This case is the first one with zero variation —,—,—,—,—,— t —,— *— ,—,—,—,—,— 

of G. Also, the exact values of “Hubble-like” 10 “ 

parameters (H = —2 h) in contrast to the 

previous case have rather simple and Figure 5: The variation of G parameter for 

compact form. m — 11. 


3. for m = 12: 


Vax(l) 


the variation of G parameter is minimal for 
1 = 11, see Figure 6, 

H « 0.2264080186, 

h ss -0.1852491999, 

Int w -0.000069032, 

Var » 0.000304899. 


Here all our four conditions are satisfied. 
The variation of G parameter is non-zero 
and the volume of the internal space is 
decreasing. 



n 1 1 1 r 

14 16 


"1 

18 


Figure 6: The variation of G parameter for 
m = 12 . 
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4. for m = 13: 


the variation of G parameter is minimal for 
1 = 9 and the “internal space” parameter is 
positive, see Figure 7, 


H « 0.2039802, h « -0.2261006, 

Int « 0.0048967, Var « 0.0240058, f °' 2 ~ 

> 

0,15 - 

for 1 = 8 the variation of G parameter is 
slightly higher, but the “internal space” pa¬ 
rameter is negative, see Figure 8, 


H« 0.1942063, h« -0.2498379, ^ I 1 T T 1 T 

1 

Int « -0.0052611, Var « 0.0263919. 

Both cases are excluded by G-dot Figure 7: The variation of G parameter for 

restrictions. m — 13 . 


5. for m = 14: 


Var(l) 

0,5 —1 

the variation of G parameter is minimal for 
1 = 7, see Figure 8, 

H « 0.1822582965, 

h w -0.2863787788, * : 

0 , 2 - 

Int « 0.000189810, 

Var » 0.00104143. °' 1 " 




A 


for m=14 


The variation of G parameter exceeds our °4— 1 — 1 — 1 —*—>—,—,—,—,—,—,—, 

limits, and the condition of the volume con- 5 1 10 

traction of the “inner space” is not met. 

Figure 8: The variation of G parameter for 
m = 14. 
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6. for m = 15: 


Var(l) 

the variation of G and the “internal space” 0 ,s- ♦ 

parameters are zero for 1 = 6, see Figure 9, 



Int = 0, 


Var = 0. 


0,2 - 


for m=15 


This is the second case with zero variation of 
G. The exact values of “Hubble-like” param¬ 
eters (H = —\h) have simple and compact 
form. 


0 


Figure 9: 
m = 15. 


1— t- 

5 


1 


"T 

10 


The variation of G parameter for 


Now we will reverse our method and look for the solutions with minimal variation of 
G for fixed l instead of m. The calculations lead to the following rule: the lesser is l the 
greater is an appropriate m which gives the minimum of Var . As we consider l > 3 and 
for l = 6 the solution with minimal variation of G parameter is already found, we should 
examine only three cases: 

1. for 1 = 5: 


Var(m) 

o.s n • 

the variation of G parameter is minimal for 0 , 4 - • • 

m = 17, see Figure 10, 

H « 0.1447364880, j°' 3 ~ 

h ps -0.4041874693, 

Int «0.005373486, * • 

0 , 1 - 

Var » 0.017555134. • . 

None of the desirable conditions are satis- a 2 o 

fied. 

Figure 10: The variation of G parameter for 
l = 5. 
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2. for 1 = 4 : 


Var(m) 

0,5 


- the variation of G parameter is minimal 
for m = 20, see Figure 11, 


- H « 0.1220672556 

- h « -0.5176845111 

- Int « 0.004405301 

- Var « 0.03608913 

The amount of variation is too high and the 
“internal space” parameter is positive. 


v 


V 


V 


20 25 

m 


Figure 11: Variation of G for l — 4. 


3. for 1 = 3: 

the variation of G parameter is minimal for 
m = 28, see Figure 12, 

H « 0.09202826388, 

h « -0.7765606872, 

Int « 0.000272078, 

Var » 0.00295646. 

The “internal space” parameter is also pos¬ 
itive and the variation of G parameter ex¬ 
ceeds the limits imposed. 


Var(m) 

0,35-, 


0,3- 

0,25 - X 



0 , 1 - 

0,05 - 


25 30 


Figure 12: The variation of G parameter for 
l = 3. 


Thus, in this subsection we have obtained cosmological solutions for m > 9, which 
satisfy all four conditions for the following cases: 


1. m = 11,1 = 16 (zero variation of G); 

2 . m = 12,1 = 11; 

3. m = 15,1 = 6 (zero variation of G). 
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It should be noted that for m — 3 and l = 2 the solution with H « 0, 750173 and 
h ~ —0, 541715 was found earlier in [16]. For this solution we have a contracting “internal 
space” bnt the variation of G is a huge one (G/G is of Hubble parameter order). Recently, 
an exact analytic form of this solution was obtained in ra- 

4 Conclusions 

We have considered the (n + l)-dimensional Einstein-Gauss-Bonnet (EGB) model. By 
using the ansatz with diagonal cosmological type metrics, we have found solutions with 
exponential dependence of scale factors with respect to “synchronous-like” variable r. 

In cosmological case (w = —1) these solutions describe an exponential expansion 
of “our” 3-dimensional factor-space with the Hubble parameter H > 0 and obey the 
observational constraints on the temporal variation of the effective gravitational constant 
G. Any solution describes (m — 3 + l) -dimensional “internal space”, which is anisotropic: 
it is expanding in (m — 3) dimensions with the Hubble rate H > 0 and contracting in l 
dimensions. 

These solutions were found (in numerical or analytical forms) for the following cases: 

1. m = 9,1 > 2680 (variation of G tends to 0 as l — > +oo); 

2. m = 11,1 = 16 (variation of G is zero); 

3. m = 12,1 = 11; 

4. m = 15,1 = 6 (variation of G is zero). 

Thus, we have shown that it is possible in the framework of EGB model to describe 
the accelerated expansion of the three-dimensional factor space with sufficiently small (or 
even zero) value of the variation of the effective gravitational constant G. For the case 
w = 1 we have obtained by product a family of static configurations which may be of 
interset within some other possible applications. 

Here we have considered a gravitational model in more than 4-dimensions. In such 
a case the Gauss-Bonnet term gives non-trivial contributions to the generalized Einstein 
held equations. In particular we have shown that there are cosmological solutions in agree¬ 
ment with observations when “projected on the (3 + l)-dimensional physical space-time. 
For the sake of simplicity, we restrict ourselves to vacuum solutions in multidimensional 
gravity with the Gauss-Bonnet term . Such ansatz may be considered as a part of a 
general “geometrical program aimed at the explanation of dark energy in 4-dimensional 
space, e.g. by using extra dimensions and modihed equation of motions just without mat¬ 
ter sources. This is a first step. The inclusion of matter sources (e.g. anisotropic fluid) 
will be the next step, e.g. as a subject of a next publication. 
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